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FLOW AND PRESSURE FIELD ANALYSIS OF PARALLEL GROOVE GEOMETRY 
FOR AN INCOMPRESSIBLE UUID WITH CONVECTIVE INERTIA EFFECTS 
byJohnZuk, Lawrence P. Ludwig, and Robert L. Johnson 
Lewis Research Center 

SUMMARY 

A set of two-dimensional equations that includes both convective inertia and viscous 
forces were formulated for parallel groove geometries. The equations are solvable for 
the velocity and pressure distribution and optimum geometry by using numerical methods 
on a high-speed digital computer. 

The mathematical model consisted of a smooth flat plate moving relative to a fixed 
parallel groove -ridge surface; both the smooth plate and the grooved plate were of infi- 
nite horizontal extent with the clearance filled with an incompressible, viscous, homo- 
geneous, isothermal, Newtonian fluid. A modified Reynolds number evolved that is based 
on a groove -ridge characteristic length parallel to the smooth plate velocity vector. The 
analysis revealed that convective inertia can only be neglected if the modified Reynolds 
number is much less than 1 (creeping flow regime). 

Application of the modified Reynolds number to four published experimental results 
reveals that a significant improvement in bearing number or sealing coefficient occurred 
when the modified Reynolds number was near 1. The detectable improvement in the 
bearing number or sealing coefficient had previously been attributed to the onset of tur- 
bulence; therefore, it is argued that the strong influence of convective inertia, and not 
turbulence, is the main cause of this improvement. 

The consideration of convective inertia effects reveals that fluid density p, aspect 
ratio b/h Q , and number of groove -ridge pairs N are additional parameters of pressure 
generation dependence that are not implied in the creeping flow analysis. 

The analysis presented herein provides a unified approach to parallel groove opera- 
tion, which was previously covered by the creeping flow analysis and a semiempirical 
turbulent flow analysis. No restriction was placed on the net flow situation between the 
smooth and grooved plates. Therefore, the analysis and resulting equations apply equally 
to hydrodynamic seals, bearings, and pumps. An arbitrary groove -ridge shape is also 
permissible in solving the equations. Consideration of convective inertia effects explains 
the principle of viscoseal axial pressure generation and attainment of zero leakage. 



INTRODUCTION 


The parallel groove geometry has many potential applications for use as a pressure 
generation device. For example, the parallel groove has been applied in spiral groove 
thrust bearings (ref. 1), spherical and conical spiral groove bearings (ref. 1), helical 
groove pumps (ref. 2), and helical groove seals (viscoseal, see ref. 3). 

The latter device (viscoseal) is potentially useful for sealing liquid metals, such as 
mercury, potassium, and sodium, in space electric power generation systems that re- 
quire 1 to 3 years of unattended operation and near zero leakage (refs. 4 and 5). As 
pointed out in references 5 and 6, the viscoseal has inherent reliability and long life 
because of the absence of solid surfaces in rubbing contact. 

The analysis of the spiral groove bearing in reference 1 indicates a higher load carry- 
ing capacity than the step or tilting pad bearings. An analysis of the helically grooved 
herringbone journal bearing is found in reference 7 for a perfect gas including the effects 
of small eccentricities. However, an infinite number of groove -ridge pairs was assumed 
for a finite -diameter journal bearing model in order to make the pressure gradients 
(which are discontinuous at the groove -ridge interface) a smooth varying function. 

To date, analyses of parallel groove geometry have been limited to creeping flow 
solutions where convective inertia effects have been neglected. Convective inertia is the 
fluid mass acceleration due to a spatial velocity change, for example, a change in ve- 
locity due to an obstacle in the flow path. This is contrasted to an inertia effect caused 
by a temporal velocity change, for example, a variation in velocity at a point in a bearing 
due to shaft runout. The creeping flow analysis is restricted to flows where the modified 
Reynolds number is much less than 1; that is, 

RS = Rej/^-f « 1 


or 


Re < 0. 01 


where 
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Creeping flow analyses for pumps, screw extruders, seals, and thrust and journal 
bearings have been published by many authors (e. g. , refs. 1 to 26). In appendix B is 
given an outline and discussion of the published creeping flow solutions. The main 
criticism of the creeping flow solution (generally used in the laminar flow regime) is that 
its validity is restricted to Reynolds numbers below most engineering applications and 
that zero net flow (seals) enters the analysis as an a priori assumption. 

A number of semiempirical formulations have been devised (refs. 4, 18, 27, and 28) 
for the operational region beyond the creeping flow regime (turbulent regime), but they 
provide no insight on the effect of various parameters and give no formal approach to 
optimization of the geometry. The main criticism of these semiempirical analyses is 
that the existance of the turbulent regime is assumed; in short, turbulence is not neces - 
sary to explain reported experimental results and may not even exist. Even if turbulence 
does exist, convective inertia effects are thought to be far more significant. 

Kettleborough (ref. 29) numerically analyzed the slider bearing with inertia, turbulent, 
and viscous terms considered. When inertia only was considered, the results were in 
qualitative agreement with published (turbulent-attributed) slider bearing experimental 
results. Kettleborough concluded that the turbulence term did not greatly affect the 
operation of the slider bearing. 

Golubiev (ref. 24) applied centrifugal pump similarity laws to the helical groove seal. 
This case considers an inviscid fluid whereby the pressure generation is found from 
Bernoulli's equation where AP varies as the square of the velocity. Due to the flow 
being in a confined space, and since the viscous force contributes to maintaining the gen- 
erated pressure gradient, the viscous force cannot be neglected. 

Unfortunately, the creeping flow optimum geometry has been extended, almost 
universally, to be valid for operation beyond the creeping flow regime. Since the convec- 
tive inertia force is the main driving force, a complete reevaluation of the optimum 
geometry including the groove -ridge shape must be investigated before conclusions on 
viscoseal performance can be made. 

The objective of this report is to examine the flow field and pressure generation in 
parallel groove geometries by including convective inertia effects. Since the groove steps 
cause convective acceleration effects of the fluid, the edges of the grooves are not ne- 
glected. A flat -plate model of infinite horizontal extent will be analyzed. This is a rea- 
sonable physical approximation because the film thicknesses (clearance gap between two 
relatively moving surfaces) are usually very small in comparison to the radius of curva- 
ture and the circumferential length of the pump, bearing, and seal rotors (screw ex- 
truders may be an exception to this). The significant terms for the study are found by a 
formal ordering procedure and by examining the physics of the problem. End effects are 
neglected because this is an infinite flat -plate case. (Refer to refs. 1 and 30 for creep- 
ing flow end effects and ref. 6 for experimental investigation of end effects. ) 
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No restriction is applied as to net flow within the parallel groove geometry; there- 
fore, the analysis applies to hydrodynamic seals, bearings, and pumps. 


SYMBOLS 



L* 

N 


n 


AP 

p* 

P c 

Pr 


Q 

Qf 

Q r 

R 


Re 

T> * 

Re 


ridge or land width 
groove width 

specific heat of fluid at constant pressure 

clearance between ridge and smooth flat plate, or radial clearance 
body force in x -direction 
body force in y -direction 
body force in z -direction 

dimensionless factor, function of helical geometry 
step or edge height 
thermal conductivity of fluid 

characteristic length, length of groove -ridge pair in smooth plate relative ve- 
locity vector direction 

axial length, normal to plate velocity 

number of helix starts 

integer 

pressure differential 

static pressure at reference state 

static cavity pressure 

Prandtl number 

net volume flow rate 

volume flow rate in plate relative velocity direction 
volume flow rate normal to plate relative velocity direction 
radius of rotor 
Reynolds number 

modified or reduced Reynolds number 
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T* temperature at reference state 

t time 

U smooth flat -plate velocity 

u velocity in x-direction 

V characteristic velocity in the y -direction 

V fluid velocity vector 

V pressure flow velocity in x-direction 
P> A 

V„ „ pressure flow velocity in z -direction 
v velocity in y -direction 

w velocity in z -direction 

x coordinate along ridge -groove 

y coordinate across film (between plates) 

z coordinate across ridge -groove 

a angle between relative velocity direction and parallel groove -ridges or helix 

angle 

£ film height 

7j coordinate normal to plate relative velocity vector 

A bearing number 

XG empirical sealing parameter (defined in ref. 5) 
ju absolute viscosity of fluid 

v kinematic viscosity of fluid 

v* kinematic viscosity of fluid at reference state 

| coordinate in direction of plate relative velocity vector 

p density 

p* density at reference state 

a shape parameter 

V vector differential operator Del or Nabla 
2 

V Laplacian operator 
Subscripts: 

c based on clearance or film thickness 
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g groove 

L based on characteristic length 

r ridge or land 

e plate relative velocity direction 

r] normal to plate velocity direction (axial length) 


BASIC MODEL 

The rectilinear Cartesian coordinate system was selected for mathematical tract - 
ability. A parallel groove -ridge plate geometry with a smooth flat plate moving relative 
to the parallel groove -ridge plate, both of infinite horizontal extent, was chosen as the 
basic model (see fig. 1). 

Conceptually this model can be thought of as an infinitely long smooth rotor with a very 
large diameter and a concentric grooved housing. Actually this infinite grooved plate 
model is the limiting case of a helical grooved cylinder where R — °° or c — 0. If 
c/R < 0. 01, this model should give a good qualitative picture of the flow field and pres- 
sure generation and significant trends for optimum geometry. For example, a practical 


oo 




(b) Cross-sectional view along z-axis illustrating relative 
motion of smooth flat plate with respect to fixed parallel 
groove-ridge plate. 

Figure 1. - Mathematical model of parallel groove-ridge geometry of infinite horizontal extent. 


6 



End effect - 


Cavity 

pressure 


Homogeneous 

fluid 


Ilk 

§4 

tf'Mk 

f§| Reg ion 

^analysis 

Invalidity 


Liquid-gas 

mixture 


■*— Scavenging or 
nonwetted region 


Average pres- 
sure ridge 
leading edge 
Average pres- 
sure ridge 
trailing edge 


-Ambient pressure 


Axial distance 




High pres- 
sure end 


WL 




■? ‘ ■ ’-3 



Pressure 

tad ' V ; *v 


C-66-2693 

Figure 2. - Region of analysis validity (homogeneous fluid) of a grooved housing viscoseal. Reynolds - 
number based on clearance or film thickness, 1700; modified or reduced Reynolds number, 19; 
400-microinch shaft rotational movement during photographic exposure; shaft speed, 8000 rpm; 
sealed fluid, water. 



High pres- 
sure end 


I Discontinuous fluid film J 

\ (mixed air and liquid) | 

Figure 3. - Viscoseal with grooved housing operating in a regime of water-air mix- 
ture (gas ingestion). Analyses does not apply for this discontinuous fluid film 
case; 230-microinch shaft rotational movement during photographic exposure; 
shaft speed, 5000 rpm ; sealed fluid, water. 

geometry is a 2 -inch -diameter rotor with a 0. 005-inch concentric clearance, 
c/R = 0.005. 

The clearance between the infinite plates is filled with any incompressible fluid, and 
it should be noted that for this infinite horizontal extent case there are no end effects 
(in practical engineering design, of course, end effects must be accounted for). Since 
experiments at Lewis show that the end effect zone decreases with increasing rotor speed 
and also that it is less than a groove width, this analysis should be valid in the bulk of 
the homogeneous fluid (see figs. 2 and 3). 

In the flat -plate case, the boundary conditions on a groove -ridge plate in motion are 
constant; however, it will be advantageous to have the smooth flat plate move relative to 
the grooved flat plate. This will eliminate the temporal inertia terms and reduce the 
subsequent independent variables by one. 

A word of caution for extension of the results to a case where the groove -ridge pairs 
on a finite shaft are rotating relative to a fixed smooth housing. The boundary conditions 
are different when the smooth housing is stationary and the grooved rotor is moving, and 
when the housing is grooved and the rotating shaft is smooth. When the rotor is grooved, 
the surface speed is different at each point along the step or edge of the moving boundary. 
Thus this analysis should be more valid for a smooth rotor and grooved fixed housing 
(from ref. 5 the grooved housing appears to be the desired mode of operation) with a 
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large radius of curvature and/or small clearance. 

When referring to figure 1, it is noted that by the choice of coordinate system the 
groove -ridge pairs appear only in the y, z -plane. Later it will be shown that this choice 
of coordinate system will be necessary to reduce the analysis to a two-dimensional 
problem for the flow field. It should also be noted that the drag force has been resolved 
into its components in this coordinate system and there is no equilibrium of forces in the 
z -direction with inertia (nonlinear) effects neglected. 


BASIC EQUATIONS 


The Navier -Stokes equations for a homogeneous, incompressible, laminar, 
Newtonian fluid are the following (see ref. 31): 
x -Direction: 


y -Direction: 


z -Direction: 


Du 3P /a 2 u 3 2 u 0 2 u) 


Dv _ 3P /0 2 v 9 2 v 0 2 v’ 

P — = F - — + Ml + + 

Dt y 0y L„2 -,„2 -,_2 


\3x 0y 


0Z 


P Dw = 0P /a?w 0?w d_K\ 
Dt Z 9z \0x 2 0y 2 0z 2 


( 1 ) 


( 2 ) 


( 3 ) 


The Eulerian or spatial derivative is D/Dt, which is composed of the local acceleration 
(temporal velocity change) and the convective acceleration (spatial velocity change): 


D 

Dt 


0 0 0 0 
— +u — + v — + w — 

0t ) y 0X 0y 0Z, 

S y y~ ' 

Local, unsteady, or Convective 

nonstationary accel- acceleration 

eration 
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The incompressible continuity equation is 


0U 0V 0W 

— + — + — = o 

0X 0y 02 


The assumptions for this analysis are as follows: 

(1) For steady flow, 


( 4 ) 


3u 

0V 

0W 

0t 

“ 0t ” 

0t 


(2) For no body forces, 



This means there are no electromagnetic fields present and gravitational effects are 
negligible. 

(3) Isothermal and equal temperature plate surface conditions are assumed. 

(4) Viscous (frictional) heating is negligible. Therefore, the thermophysical prop- 
erties are constant. This appears to be a good assumption for order ten (water) and 
smaller (liquid metals) Prandtl number Pr fluids: 



In experimental setups there is an inherent unsteadiness due to vibrations, shaft 
runout, and eccentricity of rotor with respect to housing and misalinement. These fac- 
tors will influence the degree of correlation between the analytical and experimental re- 
sults. 

The conservation of momentum and mass equations were nondimensionalized and a 
formal ordering procedure was used to determine the relative magnitude of the terms. 
Consequently, in a formal way the negligible terms were found (see appendix A). 

An important parameter, the modified or reduced Reynolds number, was found: 


Re* = Re L 



2 


= Re 


c 



if c/L < 1 
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This modified Reynolds number gives the relative magnitude of the inertia forces to the 
viscous forces; that is, 


_ * Inertia forces 

Re = 

Viscous forces 

Thus, inertia effects can be neglected only if Re* £ 0. 01 

After the formal ordering procedure (see appendix A), the three-dimensional flow 
field equations are as follows: 

Conservation of momentum: 


9u 9u 9u 1 9P 9 2 u 

u — + v — + w — = + V — — 

9x 9y 9z p 9x ^2 


(x -direction) (6) 


9P 

3y 


0 


(y -direction) 


(7) 


o 

9w 9w 9w 1 9P 9 w 

u — + v — + w — = + v 

9x 9y 9z p 9z gy2 


(z -direction) 


( 8 ) 


Conservation of mass: 


9u 9v 9w 
— + — + — 
9x 9y 9z 


= 0 


(4) 


Fully developed flow is assumed in the x-direction. This means that V and 9P/9x 
are not functions of x and that the pressure gradient in the x-direction (parallel to the 
groove) is a constant (verified experimentally in ref. 5). In figure 2, it is observed that 
this assumption is valid only where the pressure profile at the ridge leading or trailing 
edges in the axial direction is linear. Gas ingestion results in a nonlinear profile (see 
figs. 2 and 3). 

With the assumption of fully developed flow along the groove- ridge direction (x- 
direction), the flow field equations become two -dimensional: 

9u 9u 1 9P 9 2 u . 9P . . , Q v 

v — + w — = + v where — = constant (9) 

9y 9z p 9x g 2 9x 
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( 10 ) 


0P 

3y 


o 


0W dw 1 dP d^W 

V + W = + V 

3y 3z p 9z 9 2 


Conservation of mass: 


3v + dw _ o 
0y 0z 


The boundary conditions for a rectangular groove are as follows: 

( 1 ) 


v — 0 

u = U cos a 


> at y = 0 


(ID 


( 12 ) 


w = -U sin a J 


( 2 ) 


u = v = w = 0 at c^y^c + h Q (on side of groove) 
u = v = w = 0 at y=c (on ridge) 

u = v = w = 0 at y=c+h Q (on groove root) 

The boundary conditions can be generalized for an arbitrary shaped groove -ridge 
geometry (see fig. 4): 

( 1 ) 


v — 0 


u = U cos a 


> at y = 0 


w = -U sin a 


( 2 ) 


u = v = w = 0 at y = h(z) 
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Figure 4. - Cross-sectional view along z-axis 
illustrating arbitrary groove- ridge pair shape. 


Equations (9) to (12) can be reduced in terms of 
the dependent variables u, v, and w by vanishing 
the pressure terms. The pressure terms vanish when 
differentiating equations (9) and (11) with respect to 
y and when noting that 3P/0y = 0 (eq. (10)). In the 
process the order of the equations is increased by 
one. The x-direction momentum equation becomes 


0V 0U 0W 0U 0^U 0^U 0^U 

+ + w + v = v 

0y 3y 3y sz dydz 9y 2 dy 3 


The z -direction momentum equation becomes 


2 2 3 

0V 0W 0W 0W 0 W 0 W 3 w 

+ — — + w + v = v 

0y 0y 0y 0z 0y0z .2 3 

oy a y 


Since the continuity equation remains the same, 


0V 0W . 
— + — = 0 
0y 0z 


then the z -direction momentum equation becomes 


w 


.2 a 2 

0 W 0 W 

+ V 

0y0z ay 2 



(13) 


(14) 


( 12 ) 


(15) 


These equations (12), (13), and (15) can be solved for the velocity components. Once 
the velocity components have been obtained, the pressure gradients can be found from 
equations (9) and (11). Then the axial pressure gradient 0P/0L’ can be found from 


0P 0P . / 0P\ 

= — sm a + — cos a 

0l* 0 X W r+g 


(16) 


where L’ is the axial length (see fig. 5). Pressure distribution in the x- and z-directions 
is shown in figures 6 and 7. Therefore, by examining the analysis and the derived equa- 
tions, the axial pressure gradient 0P/0L' is a function of (a, U, p, v), where a is a 
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(b) Across ridges and grooves (no restric- 
tion that profile must be linear). 

Figure 7. - Pressure profiles (qualitative graphical representation from experimental data). 
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shape parameter that is a function of (b/h Q , c/h Q , a/b, a). 

Since convective inertia effects are important, new geometric parameters appear 
which are not implied in creeping flow theory. Most notable of these is the aspect ratio 
b/h Q , which implies the degree of convective effect. When the b/h Q and a/b ratios are 
known, N, the number of helix or groove starts, can be found. Thus the number of 
groove or helix starts also is an important parameter. Since inertia is important, a 
high density is desired for pressure generation. 


DISCUSSION 

Discussion of Equations 

This analysis is valid only for a homogeneous fluid. Figures 2 and 3 clarify the 
region of validity. In some parallel groove geometries (e. g. , the viscoseal) where there 
is a liquid-gas interface, there is a problem with gas ingestion (ref. 5). Gas ingestion, 
which is a condition of liquid and gas mixture throughout the seal, is shown in figure 3. 
The analysis is not valid for this condition because the fluid is nonhomogeneous. Further, 
the analysis is not valid for a scavenging area and a region consisting of a mixture of gas 
and liquid (see fig. 2). 

It should also be pointed out that the analysis is not valid in the end effect region at 
the high pressure end of the seal (see fig. 2) where the groove connects with the cavity. 
Since the appropriate equations are nonlinear, the pressure is not a harmonic function; 
thus, the usual analog methods of finding end effect are not available as found in creeping 
flow regions (see appendix B). 

The equations are unique in many ways. Mathematically, the problem has been re- 
duced to two dimensions (y- and z -directions) but all three velocity components are 
present. Also, there is a pressure gradient in the third dimension (x -direction), but 
fortunately this gradient is a constant. 

The velocity component v is no longer negligible as it was in the creeping flow solu- 
tion. Physically, v cannot be negligible because of the presence of the edge effect which 
will result in a component of velocity in the y-direction. 

Since there is no restriction on net volume flow, the pressure generation prediction 
equations should equally apply to hydrodynamic seals, bearings, and pumps. 

When further examining the physics of the problem, it is seen that across the 
groove -ridge pairs the unbalance of the convective acceleration force with the drag force 
results in a pressure gradient in the z -direction. Since the step or edge causes the con- 
vective acceleration, its effects are not negligible unless there is a very large aspect 
ratio and the mathematical model accounts for the edge effect. Along the groove and 
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ridge (x -direction) the drag force is in equilibrium with the induced pressure and con- 
vective inertia force. The drag force not only causes fluid motion but also aids in main- 
taining equilibrium with the pressure along the groove or ridge. 

Both equations of motion illustrate the coupling of u, v, and w in the convective 
inertia terms. This means that physically a condition such as no end leakage is possible. 
In creeping flow there is no coupling between the velocity components, which illustrates 
the shortcoming of that analysis (i. e. , end leakage is not zero). 

The assumption that the inertia terms are of the same order as the viscous terms re- 
sults in equations of the boundary layer type; however, the boundary conditions are dif- 
ferent. 

The equations are of a form that might tempt one to reduce the number of independent 
variables by one by trying to find the admissable flow field which would result in a 
similarity transform solution. This technique fails because of the characteristic lengths 
that are present in this type of confined flow. 

The equations are solvable by using high-speed digital computer numerical methods. 

A mathematical function to describe the groove -ridge shape would be too complex for a 
closed -form analytical solution. It should be interesting to find out if the step is still 
the optimum shape with convective inertia effects included. For creeping flow, Lord 
Rayleigh found in 1912 that the step was the optimum shape by using calculus of varia- 
tions. 

The equations and ordering procedure should be reexamined carefully before extend- 
ing these results beyond the assumptions used in this analysis. 


Discussion of Modified Reynolds Number 

A check was made on the significance of the modified or reduced Reynolds number 
by calculating Re* from the available published data at the claimed onset of turbulence 
(break point). In table I it can be noted that the departure in all cases from a constant 
sealing coefficient occurred when Re* was in the range of 0. 26 to 2. 2, which suggests 
that the so-called beginning of turbulent operation is really the beginning of significant 
convective inertia effects. This would be in agreement with Kettleborough (ref. 29), who 
found qualitative agreement between his inertia solution and experimental turbulent re- 
sults. 

For seals 5, 6, and 7 (see ref. 27 and table I), only the groove to ridge ratio has 
changed. As the ratio increases, the breakpoint decreases. This is reasonable since a 
larger ridge area results in a less influential convective inertia effect and delays the 
beginning of pronounced convective inertia effects. It can also be seen in table I that the 
calculated Re* are within an order of magnitude of one another for all of the experi- 
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TABLE I. - CRITICAL MODIFIED REYNOLDS NUMBER FOR ONSET OF TURBULENCE OF VARIOUS EXPERIMENTERS 


Experimenter 

Reference 

Seal 

Diameter, 

in. 

Radial 

clearance, 

c, 

in. 

Ratio of 
radial 

clearance to 
rotor radius, 
c/R 

Ridge 

width, 

a, 

in. 

Groove 

width, 

b, 

in. 

Helix 

angle, 

deg 

Characteristic 

length, 

L, 

in. 

Clearance 

Reynolds 

number, 

Re c 

Modified 

Reynolds 

number, 

Re* 

■ 

Stair 

27 

a l 

1.2430 

0.0042 

0.0068 

0.0934 

0.1596 

14.5 

1.01 

400 

1.7 




2 

1.2465 

.00235 

.0038 

.1176 

.0494 

9.67 

.98 

(b) 

(b) 




2B 

1.2420 

.0047 

.0076 

.1176 

.0494 

9.67 

.98 

(b) 

(b) 




3 

1.2461 

. 00265 

.0043 

.0828 

.0842 

9.67 

.98 

350 

.95 

Z 



3B 

1.2408 

.0053 

.0085 

.0828 

.0842 

9.67 

.98 

(b) 

(b) 

= 



4 

1.2461 

. 00265 

.0043 

.0514 

.1156 

9.67 

.98 

200 

.95 




4B 

1.2408 

.0053 

.0085 

.0514 

.1156 

9.67 

.98 

400 

2.2 




5 

1.2455 

. 00295 

.0047 

.1585 

.0607 

5. 81 

1.83 

600 

.97 

■ 



6 

1.2455 

. 00295 

.0047 

.1070 

.1083 

5. 81 

1.83 

200 

.32 

■ 



7 

1.2455 

. 00295 

.0047 

.0609 

.1600 

5. 81 

1.83 

160 

.26 

: 

King 

4 and 27 

3C 

2 

0. 0029 

0.0029 

0. 1063 

0. 1063 

3.86 

3.2 

550 

0. 50 




IB 

2 

.0020 

.0020 

.1063 

.1063 

3.86 

3.2 

400 

.25 




2E 

2 

.0025 

.0020 

.1000 

.1000 

7. 26 

1.54 

650 

1.1 

i 



5D 

2 

.0015 

.0015 

.1000 

.1000 

7.26 

1.54 

650 

.63 

! 

McGrew and 

18 

SF96-5 oil 

1 

0. 00325 

0.0065 

(c) 

(c) 

14.5 

3.14 

400 

0.41 


McHugh 













Lessley 

28 

Configura- 

2 

0. 0032 

0.0032 

0.075 

0.122 

14.5 

0.79 

350 

1.4 




tion B 













Configura- 

2 

.0032 

.0032 

.083 

.136 

4 

3.13 

500 

.51 




tion C 












a Eccentricity of rotor to housing, e = 0. 1. 
p Break point not well defined. 
c Four thread starts. 



Empirical sealing coefficient, XG (log scale) Sealing parameter, XG, dimensionless 


.4 


.1 

.08 

.06 


1 

O 
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Figure 8. - Experimental sealing parameter tor seals 5, 6, and 7 from reference 27. 
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Figure 9. - Empirical sealing coefficient as function of 
clearance Reynolds number. 


menters. It is difficult, nevertheless, to 
select the breakpoint (see fig. 8). Even in 
a range of possible breakpoint numbers, 
however, the order of magnitude will re- 
main the same. Further experimental re- 
sults will be necessary, but in general, it 
can be concluded that if Re* si, a constant 
sealing coefficient or bearing number will 
not exist. 

In figure 9 it is seen that classically 
the empirical sealing parameter XG was 
thought to be a constant until some critical 
Re c and then to vary exponentially in the 
turbulent regime. The creeping flow solu- 
tion was thought to be valid until the critical 
Re c and then a semiempirical turbulent 
theory had to be used. Actually, the onset 
of significant convective inertia should be 
gradual, as shown in figure 9. Data (from 
ref. 27) for seal 7 in figure 8 tends to show 
this. The derived equations herein pertain 
to all cases under the assumptions stated 
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and will break down only under strong turbulence or an unknown instability pnenomena. 


Principle of Viscoseal Operation 

The classical explanation of viscoseal operation is now described. Since the seal is 
’’pumping'’ the sealed fluid as fast as it is leaking out, there is no net end flow (zero 
leakage). From the basic laws of fluid mechanics, however, it is seen that the visco- 
seal axial pressure gradient results from the unbalance of the convective inertia force 
across the groove -ridge pairs and the drag force. The convective inertia force is much 
greater than the opposing drag force and results in a net pressure force. Along the 
groove -ridge direction, the drag force maintains equilibrium with the pressure force and 
convective inertia force. All three velocity components are "coupled”, and thus zero 
net leakage is achieved. This is an excellent example of the physical importance of a 
nonlinear effect. The inadequacy of the linear theory (creeping flow) has been previously 
discussed. 


CONCLUDING REMARKS 


A review of the literature on analysis of the parallel groove geometry revealed that 
either the convective inertia forces or the viscous forces were neglected (inviscid fluid 
assumption). This analysis included both convective inertia and viscous forces. The 
following two-dimensional set of flow field equations were formulated and are solvable 
for the velocity and pressure distribution and optimum geometry utilizing numerical 
methods on a high-speed digital computer: 

Conservation of momentum: 


9u 9u 1 9P 9 2 u 9P 

v — + w — = + v — = constant 

9y 9z p 9x j 2 9x 


(9) 


9P 

3y 


0 


( 10 ) 


9w 9w 1 9P 9 2 w 

v — + w — = + v 

9y 9z p 9z 9y 2 


(ID 
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Conservation of mass: 


9v 3w . 

— + — = 0 
3y 0z 

With boundary conditions at the moving smooth plate surface, 

u = U cos a 
v = 0 


( 12 ) 


w = -U sin a 


With boundary conditions at the fixed parallel groove geometry surface, 

u = v = w = 0 

From the development and analysis of the equations, the following results and con- 
clusions can be made: 

1. A modified or reduced Reynolds number evolved which gives the relative magni- 
tude of the convective inertia forces to the viscous forces: 



The characteristic length was selected to be the length of the groove -ridge pair in the 
direction of the smooth plate velocity vector. This is the principle length over which the 
driving force acts. The flow behavior over each groove -ridge pair is the same. Convec- 
tive inertia effects can only be neglected if Re* « 1, which will result in true creeping 
flow. 

2. A modified Reynolds number equal to 1 will predict a significant change in sealing 
coefficient or bearing number which heretofore was attributed to the onset of turbulence. 
This means that the convective inertia effects are the foremost contributing factor to the 
improvement in sealing coefficient or bearing number, not the turbulence effects. A 
calculation of the modified Reynolds number at the onset of turbulence of four different 
experimenters disclosed that the modified Reynolds number was near 1 (Re* = 0. 25 to2. 2). 

3. Several additional parameters of pressure generation dependence where found 
that are not implied in the creeping flow analysis. These additional parameters, which 
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arise from convective inertia effects, are the fluid density p, the groove width to depth 
or aspect ratio b/h Q , and the number of groove -ridge pairs N. 

4. Solving equations (9) to (12) by numerical analysis will permit the investigation of 
arbitrary groove -ridge shape, for example, sinusoidal, rectangular, triangular, etc. 

(It is expected that the optimum groove -ridge geometry will be different from that found 
from the creeping flow solution. ) Now a check can be made to see if a step is the opti- 
mum shape with convective inertia effects as it is for creeping flow. 

5. The analysis provides a unified approach to parallel groove operation previously 
covered by creeping flow (sometimes called laminar flow) analysis and a semiempirical 
turbulent flow analysis. 

6. No restriction was placed on the net flow situation between the plates. Therefore, 
the analysis can be used to analyze parallel groove geometries in hydrodynamic seals, 
bearings, and pumps. 

7. The viscoseal axial pressure gradient results from the unbalance of the convective 
inertia force across the groove -ridge pairs and the drag force. Along the grooves or 
ridges, the drag force maintains equilibrium with the pressure force and the convective 
inertia force. 

8. Zero leakage in viscoseal analysis can be attained through the consideration of 
convective acceleration terms (nonlinear effect). It is argued that the creeping flow 
analysis, which is a linear theory, cannot prove zero leakage. 

9. The format of this analysis can be used to evaluate qualitative convective inertia 
effects in any step or converging wall hydrodynamic seal, pump, or bearing (Rayleigh 
step, journal bearing, etc.). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 26, 1966, 

120-27-04-03-22. 
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APPENDIX A 


FORMAL ORDERING PROCEDURE FOR SIMPLIFYING BASIC EQUATIONS 

The Navier -Stokes equations for a homogeneous incompressible Newtonian fluid, 
neglecting body forces and for steady flow, are 

(Al) 
(A2) 


(A3) 


(A4) 

To find the important terms in the previous flow field equations, a formal ordering 
procedure will be used to determine the relative magnitudes of the terms. 

The terms are normalized by scaling them to their appropriate maximum charac- 
teristic value as follows: let 


/ 3u 3u 3u\ 

p I u — + v — + w — ) = 
\ 3x 3y 3z I 


ap /a 2 u a 2 u a 2 u\ 

— + fl ( + + 

8 * \9x 2 ay 2 


( dv dv 3v\ 3P 

p I u — + v — +w — ] = - — + ju 

\ dx dy dz I dy 


a 2 v + a 2 v + aV 
^3x 2 ay 2 az 2 y 


/ 3w 3w 3w\ 

p [ u — + v — + w — 1 = 

\ 3x 3y 3z I 


ap /a 2 w a 2 w a 2 wN 

— + p ( + + 

9z \ax 2 ay 2 3z 2 / 


The incompressible continuity equation is 


3u 3v 3w 
— + — + — = 0 
dx dy dz 


x = 


u = 


X 

v - 3L 

— 

z 

L cos a 

y - 

c 

z — 

L sin a 

u 

— V 

w = 

w 

U cos a 

v = — 

V 

U sin a 


The characteristic dimensions can be seen in figure 1. It is noted that a trigono- 
metric relation exists between the x and z characteristic lengths and that c/L < 1. 
The characteristic length was selected to be the length of the groove -ridge pair in the 
direction of the smooth plate velocity vector. This is the principle length over which 
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the driving force acts. The flow behavior over each groove -ridge pair is the same. 

Let p* and v* be the fluid density and viscosity corresponding to a reference state 
(P*, T*). Therefore, 


v = 


V 

v* 


1 


and 


P = 



1 


for incompressible fluids with constant thermophysical properties. Now to find the rel- 
ative magnitude of V so that all terms of the incompressible continuity equation will be 
of the same order, the incompressible continuity equation is nondimensionalized: 


9u V L 9v [ 9w 
9x U c 9y 9z 


(A5) 


To have all terms of the equation the same order of magnitude, — — must be of the order 

off. Therefore, V = U— , and thus v = — 

L U c 

Next the equations of motion are nondimensionalized. The pressure term must first 
be nondimensionalized; even though the pressure is not the driving force, it is the func- 
tion of primary importance. Thus the characteristic pressure is not extremely impor- 
tant. Since the pressure gradient is related to the rotor speed, it appears to be logical 
to nondimensionalize the pressure with respect to the dynamic pressure pU (also called 
the dynamic head): 


P 


P 



o 

It can be seen that the characteristic dynamic pressure pU has reduced the num- 
ber of nondimensional groups by one (i. e. , if P* (a reference pressure) had been used, 
another nondimensional group would be present in the equations). Also, since the incom- 
pressible case is being analyzed, the pressure appears in the dynamic role only, not in 
the thermodynamic role. 
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When the x-direction momentum equation is nondimensionalized, equation (1) re- 
sults in 


u — + v — + w 
dx dy 


dn 

1 0P V* 

" 1 >ts +C- 

L\ 2 9 2 u + 1 9 2 u 

dz 

2 dx UL 
cos a 

_COS 2 Q! 3 X 2 

z ' Sy 2 sin 2 ci! 0z^_ 


Next a Reynolds number is defined based on the characteristic length L: 


Ro L - 


UL 

v* 


where 


l> = V*V — u* 


and thus 


_ 9u _ 3u _ 
u — + v — + w 
3x dy 


011 _ 

1 

3P 1 

1 

CD 

CO 

£1 


1 3 2 u 

0Z 

cos 2 at 

3x Re^ 

_cos 2 a 


J «* + 

sin 2 a Sz 2 - 


The convective acceleration and pressure terms are of unit order; thus, the viscous 

term must be of unit order also. The largest term of the viscous forces is 

2 2 2 

— — > since f — \ <1; therefore, the (— ) should be of unit order. The mod - 

a -2 [lJ Re L \cj 


ified Reynolds number is defined by Re* = Re^ . This is the identical "reduced" 
or "effective" Reynolds number that is used in theoearing theory to show the true rel- 
ative magnitude of the inertia and viscous forces; therefore, 


_ * Inertia forces 

Re = 

Viscous forces 

(e.g. , see ref. 31). Now for curved geometries, a Reynolds number based on clearance 
has appeared. This is especially true in viscoseal analysis where 


Re 


c 


Uc 

v 
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Writing the previous x-direction momentum equation in terms of Re gives 

V 


_8u _3u^_9u 1 0P v* (c 1 0 2 u L 0 2 u c 1 0 2 u\ 

u— +v— + w — = — + — [ + + ) 

9x 0 y 9z cos 2 a 0x Uc \ L cos 2 * 0X 2 c df L sin 2 a 0z 2 / 


or 


— 0u — 0U — 0U 

U — + V — + W — : 

0x 0y 0z 


cos 2 a 


0P ( 1 

0x Re, 


cos 2 * 0X 2 


0 2 u | L 0 2 u 


c 1 d u 

c rs—2 L . 2 2 

3y sin a? 9z , 


O o 

Again it is seen that the d v/dy term is the largest. 
Thus l/Re c (L/c) should be of unit order. Now 


Re* 



which leads to 


Re* = Re^ 




It does not matter which Reynolds number is used to find the relative inertia force to 
viscous force ratio provided that c/L < 1, which corresponds to the conditions 


and 


? 2 a 2 
— » — 

3y 2 3z^ 



The proper x-direction momentum equation then becomes 
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(A6) 


3u 3u 3u 1 3P 3^u 

u — + v — + w — = + p. 

3x 3y 3z p 3x 9 2 


It is assumed that a does not go to 0° or 90° for the formal ordering procedure to be 
valid. 


Example 

A numerical example will be used to illustrate the validity of neglecting the 3^/9z 2 
2/2 

and 3 /3x viscous terms. A typical viacoseal has the following construction and 
operating conditions: smooth 2 -inch -diameter rotor, 0. 0035-inch concentric radial 
clearance between the ridge (housing) and rotor surfaces, 1000 rpm, and sealed fluid 
(water) at 60° F, where 

v = 17. 5X10“^ sq in. /sec 

c = 0. 0035 in. 

U = 2n-R(1000 rpm)/60 = 100 in. /sec (plate speed) 

c/R = 0. 0035 (therefore, curvature effects should be slight) 

L = Ridge width + Groove width (in circumferential direction) 

= 0.2 in. +0.2 in. = 0. 4 in. 

a = 15° 


Re L = UL/u = 23 000 
Re c = TJc/i> = 200 

Using the Reynolds number based on the characteristic length 



1 3 2 u 

2 TJ + 

cos a 3x 




1 3^u 

9 _9 

sin a 3z 



Using the Reynolds number based on the clearance c gives 


L 

13 


yields 

a 2— a 2-\ 

000 ^+ 15 ^] 
Sy 2 3Z 2 / 
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1 8 2 u 

cos 2 a Sx 2 


+ 


(h\& n 

z\ i a 2 u 


u . 2 --2 
^ sm a 3z 



— + 114 — + 0.13 
83E 2 ay 2 



o p 

In both equations it is seen that 3 u/3y is the largest term. Now the modified or re- 
duced Reynolds number is 


Re* 




1.75 


Nondimensionalization of y- and z-Direction Momentum Equations 

Nondimensionalizing the y-direction momentum equation (eq. (2)) yields 


_ 3v _3v _ dv 

u — + v — + w — = 
0x dy dz 



ap l 

+ 

dy Re^c 


l a 2 v 

2 —2 + 
cos* a air 



i a 2 v 

+ 

sin 2 a 3z^ 


Again using the definition of Re* and comparing the iargest terms give 


1 


Re* 



0 ( 1 ) 


and 

/l\ 2 ap l /l\ 2 aV 
W 3y Re L \c/ 0 y2 

Since 


1 


Re 


L 


« 1 


therefore 
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££ 0 

3y 


(A7) 


This result is expected for a flat plate where there is no mechanism to generate a radial 
pressure gradient (e. g. , a centrifugal force). For small curvatures, however, the 
radial pressure gradient 3P/3y should also be a small quantity and have little physical 
influence on the axial pressure gradient. 

In a similar manner the z -direction momentum equation is nondimensionalized and 
the result is 


9w 3w 3w 1 3P 3^w 

u — + v — + w — = + M 

3x 0y 0z p 0z gy2 


(A8) 


As expected, due to dimensions of tne same relative size, the z-direction momentum equa- 
tion is the same as the x-direction momentum equation (eq. (6)), but the boundary condi- 
tions are different. 



APPENDIX B 


SUMMARY OF CREEPING FLOW ANALYSIS 

The optimization from creeping flow analysis has been formulated by considering a 
flat -plate model (unwrapped cylinders, see fig. 10 ) and the following: 

(1) General incompressible plane Couette flow is assumed in the £ -direction 
(Poiseuille flow and simple Couette flow). The £ -direction momentum equation is 

_1 9P 

0^2 p 0 | 


with the boundary conditions 


u £ = U at £ = 0 
U £ = 0 at £ = h 


where 


h = c at ridge surface 
h = c + h Q at groove root surface 
The well known solution (e. g. , see ref. 31) is 



Figure 10. - Model used in creeping flow solution. 


« =-U 
6 h 


A . A P 

2 p d£ l "hi 


(2) Poiseuille flow is assumed in the 
^-direction. The 77 -direction momentum 
equation is 

a 2 

9 w 7? _ 1 9P 
a? 2 p dri 


with the boundary conditions 
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U x - 0. 96 U 


(High pressure) 


w^ = 0 at £ = 0 



U z - 0.261) 


{Low pressure) 


Figure 11. - Resolution of plate velocity {or drag force) 
into components along and across groove or ridge. 


w = 0 at £ = h 
V 

The well known solution (see ref, 31) is 


w 


V 


= rfF) - h) 

2 pt \ 077 / 


Now the condition Q = 0 is imposed in the control volume; that is, 


Q = Q| + = fj u £ d7? d£ + jj w^d£ di; = 0 

which means that there is no net end leakage or closed channel flow. This is a restric- 
tion that limits the optimization and pressure gradient prediction to a seal only. 

Solving for pressure gradient that satisfies the Reynolds equation (ref. 32) results in 
a sealing parameter 


AP 

L* 


6juUG 


or 



6pU 



1 

A 


where A is the bearing number (ref. 32) and G is a geometry factor. See references 1 
to 26 for various expressions and values for G. 

The creeping flow analysis has some paradoxes. It can best be illustrated by re- 
ferring to figures 1 and 11 and by observing the resolution of the drag force into com- 
ponents along the groove (x -direction) and normal to the groove (z -direction). In creep- 
ing flow, the drag force is in equilibrium with the pressure force. Thus along the groove 
(x -direction) the drag force is indeed in equilibrium with the pressure force. In the 
normal groove direction (z -direction), however, the pressure and drag forces are in the 
same direction, which is also the direction of a very large pressure gradient (see figs. 

6 and 7). Therefore, an unbalance of forces exists in this direction. Physically, the 
Q = 0 restriction is not possible. Also, this point can be illustrated from £- and 
7j -direction momentum equations, respectively: 
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Notice that there is no velocity ’’coupling" between the velocity components; that is, 
u £ and are independent of one another. It is obvious that for zero leakage these 
velocity components are not independent of one another. Thus, physically as well as 
mathematically it is not possible under the restrictions of the classical creeping flow 
analysis to get zero leakage. 

Also in the creeping flow analysis the edge effects are assumed negligible. The full 
meaning of this assumption has not really been clarified; it is only close to physical 
reality when c is very small and/ or b/h Q is very large. For example, the creeping 
flow analysis is valid for spiral groove thrust bearings where c is on the order of 
500 microinches and the aspect ratio is on the order of 1000 (c/L « 1, see fig. 12). As 
b — h Q , the edge effect becomes more pronounced as the convective force becomes more 
and more important. To date only reference 30 has explicitly stated this restriction. 

To resolve the paradox of having the drag and pressure forces unbalanced in the 
z -direction, a convective inertia force is in equilibrium with the pressure and drag 
forces (a nonlinear effect). The step is causing a convective change that results in gen- 
eration of an axial pressure gradient. The role of the step is now described. 

If a concentric rotor is rotating about a stationary smooth sleeve, no axial pressure 
gradient is generated. This is the classic axisymmetric Couette flow (see ref. 33). 
Furthermore, if an axial pressure gradient is imposed (e. g. , by having the rotor trans- 
late in the axial direction) an "uncoupled” helical flow will result. Subsequently, of 
course, there will be net flow out. 

In the helical groove seal (viscoseal), the step or edge h Q is the pressure genera- 
tion mechanism. Without the step, that is, the limiting case when h n -*• 0 (rotor and 




h 0 

W 4 

'/ 


/ 


\m 1000 h 0 or 1000 c J 



Figure 12. - Example where the creeping flow mathematical model can be considered a valid 
physical model. Modified Reynolds number (Re*) much less than 1 since clearance c is 
much less than characteristic length L. 



housing smooth), there is no axial pressure gradient generated. Eccentricity will gen- 
erate a radial pressure gradient. Therefore, the effect of the groove wall is a convec- 
tive acceleration or inertia effect (nonlinear effect). To neglect the convective force 
terms means that the primary axial pressure generating mechanism is neglected. 

The creeping flow solution equations are linear and pressure is a potential function. 
The creeping motion or slow viscous flow momentum equation in vector notation is 

VP = pV 2 V 


The incompressible continuity equation says 

V • V = 0 

Taking the divergence of both sides of the equation yields 

V • VP = pV 2 (V • V) = 0 


or 


V 2 P = 0 

which is Laplace’s equation, since pressure is a scalar quantity. Consequently, the 
multitude of mathematical analog methods of solving linear equations and specifically 
Laplace’s equation can be used. Thus, solutions for end effect can readily be found 
(refs. 1 and 30). 

Muijderman (ref. 1) used conformal mapping to analyze the spiral grooved bearing. 
Again, this technique can be readily applied because pressure is a potential function. 
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